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Objectives 

Our  Phase  II  work  seeks  to  produce  software  for  efficient  and  accurate  evaluation  of  elec¬ 
tromagnetic  fields  on  and  around  aircraft,  with  application  to  antenna  prediction  problems. 
Our  methods  are  based  on  two  main  elements  we  introduced  recently:  a)  A  new  class  of 
rigorous  high-order  integrators  applicable  to  problems  of  high  frequency,  and  b)  High-order 
surf  ace- re  present,  at  ion  methods  based  on  a  novel  Fourier  continuation  technique. 


Figure  1:  Left:  Original  low  quality  surface  rendering.  Center:  processed  surface  by  means 
of  the  continuation  method  surface  representation  algorithm.  Right:  Solution  of  a  200 
wavelength  problem  of  scattering  by  a  surface  resulting  from  the  surface  representation 
approach. 


Status  of  Effort 


Having  tackled,  over  the  first  year  of  this  effort,  problems  involving  high-performance  par¬ 
allel  implementations,  and  development  of  a  quality  Graphical  User  Interface  (GUI),  during 
the  second  year  of  this  work  we  turned  to  the  four  remaining  development  items  listed  be¬ 
low:  over  the  forthcoming  six  months  of  the  effort  (on  a  no-cost  extension  of  the  present 
contract)  we  will  integrate  all  solvers  into  the  previously  developed  GUI,  and  we  will  ad¬ 
dress  a  number  of  issues  concerning  singular  geometries  and  curved  antenna  problems,  as 
discussed  in  detail  in  this  report.  The  four  main  items  that  were  addressed  during  this 
performance  period  are  as  follows:  1)  Solution  of  problems  of  Electromagnetic  Scattering  by 


Open  Surfaces  (we  believe  that  the  contributions  reported  below  in  these  regards  provide  for 
the  first  time ,  a  stable  and  accurate  solver  for  open  surface  scattering  in  electromagnetics 
(other  previous  open  surface  solvers  were  limited  to  cases  in  which  physical  fields  do  not 
tend  to  infinity  at  the  edge;  approximations  found  in  the  engineering  literature,  in  turn, 
do  not  use  sufficiently  refined/accurate  solvers  and  do  not  produce  reliable  approximations; 
2]  Enhancement  and  extension  of  the  previously  developed  Geometry  Representation  Soft¬ 
ware:  3)  New  and  simplified  of  methodologies  for  treatment  of  Corner-  and  Edge-scattering; 
4}  Integration  of  surface  representation  and  solver  capabilities,  with  application  to  Solution 
of  scattering  problems  for  large,  realistic  engineering  structures ,  and  5)  Curved  wire  antenna 
solvers. 


1  Electromagnetic  Scattering  by  Open  Surfaces 

The  goal  of  this  aspect  of  our  work  is  to  build  an  efficient  and  accurate  algorithm  to 
solve  electromagnetic  scattering  by  open  surfaces;  as  mentioned  above,  we  believe  that 
the  contributions  reported  below  in  these  regards  provide  for  the  first  time t  a  stable  and 
accurate  solver  for  open  surface  scattering  in  electromagnetics  (our  previous  open  surface 
solvers  wore  limited  to  the  acoustic  equations;  other  previous  contributions,  in  turn,  do  not 
use  sufficiently  refined/ accurate  solvers  to  even  approximate  the  conditioning  difficulties 
inherent  in  this  problem). 


Figure  2:  Left:  model  open  surface  scatterer:  annulus;  note  the  arrow  indicating  the  prop¬ 
agation  direction  of  the  incident  field,  and  the  radius  indicating  the  polarization,  that  is, 
the  orientation  of  the  electric  field  vector.  Right;  resulting  current  on  the  annulus.  Note 
that,  in  spite  of  the  symmetry  of  the  incidence,  the  current  is  not  circularly  symmetric — as 
a  result  of  the  electrical  polarization. 

The  key  property  of  the  boundary  integral  formulation  of  Maxwell  equations  is  that  the 
solution,  i.e.  surface  current  densities,  exhibit  two  types  of  singularity  near  the  open  edge 
of  the  boundary;  the  tangential  component  of  the  current  density  tends  to  infinity  and  the 
normal  component  decays  to  zero  as  the  edge  is  approached.  More  precisely,  the  tangential 
(Ju)  and  the  normal  ( Jv )  components  are  given  by  the  products  of  known  singularities 
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Figure  3;  Electromagnetic  scattering  by  an  annulus  with  inner  radius  0.5  and  outer  radius 
1.  Ey  is  displayed.  The  incident  plane  wave  has  direction  (0,071)?  polarization  (0,1,0). 
The  wavenumber  is  i r,  the  corresponding  wavelength  is  the  same  as  the  aperture  diameter. 
On  the  annulus,  the  tangential  component  of  the  surface  current  density  is  displayed,  which 
tends  to  infinity  near  the  open  edge. 


(1  fyfd  and  \fd)  and  functions  (au  and  av)  smooth  up  to  the  boundary; 

Ju  —  du/yfd  (1) 

tJy  =  Cly  "s/d  (2) 

where  d  is  the  distance  from  the  open  edge.  The  electric  field  integral  equation  in  terms  of 
au  and  av  becomes 


p  zfcv  p  %  $£7* 

Twa  =  k2v  X  /  y)Wa(y)dSv  +  vxVx  /  div(Wo{if))dS;  (3) 

Jy  4tt r  Jy  47rr 

-  iky  x  Einc 

where  n(:7\?/)  is  the  projection  operator  from  the  tangent  space  at  y  onto  the  tangent  space 
at  x ,  and  where  W  —  ^  j .  The  operator  T\y  consists  fundamentally  of  two  weighted 

single  layer  potentials  with  appropriate  pre-  and  post- applications  of  surface  divergence  and 
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Figure  4:  Left:  model  open  surface  scatterer:  annulus;  note  the  arrow  indicating  the  prop¬ 
agation  direction  of  the  incident  field,  and  the  radius  indicating  the  polarization,  that  is, 
the  orientation  of  the  electric  field  vector.  Right:  resulting  current  on  the  annulus.  Note 
that,  in  spite  of  the  symmetry  of  the  incidence,  the  current  is  not  circularly  symmetric — as 
a  result  of  the  electrical  polarization. 


gradient.  We  applied  typical  pseudo-spectral  differentiation  schemes  for  the  divergence  and 
the  gradient.  The  main  numerical  difficulties  arise  when  we  evaluate  the  integral  with  the 
singular  scaling  factors  W.  We  established  an  evaluation  scheme  for  the  weighted  single 
layer  potential  as  follows: 

•  Let  [—1,1]  x  [—1,1]  be  the  parameter  space  of  each  patch  of  the  surface,  and  let  the 
open  edge  correspond  to  the  lines  v  =  ±L  We  introduced  Che  by  she  v  grid  for  v.  that 
is,  we  applied  a  change  of  variable  v  =  cos (t)  and  utilized  a  uniform  grid  for  t.  Then, 
</d=y/l-  y2  —  sin(£).  and  the  scaled  density  and  its  divergence  become 


(4) 


d\v{Wa)  = 


1 

sin(£) 


(du  {au)  +  dt(av)  sin(£)  +  av  cos (£)) 


(5) 


The  common  1/  sin(£)  term  is  canceled  out  by  the  sin(t)  term  arises  from  the  Jacobian 
of  the  change  of  variable.  Thus,  the  integrands  become  smooth  and  even  functions  of 
t,  except  for  the  1  jr  singularity  of  the  kernel, 

•  The  kernel  singularity  can  be  handled  by  changing  the  integral  in  polar  coordinates  as 
our  previous  scheme  for  regular  closed  surface  cases.  However,  unlike  the  closed  surface 
cases,  the  polar  integral  for  open  surfaces  involves  two  serious  numerical  difficulties: 
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Figure  5:  Side  view.  Left:  Original  low  quality  surface  rendering.  Right:  processed  surface 
by  means  of  the  continuation  method  surface  representation  algorithm. 


-  For  closed  surfaces,  the  integrand  vanishes  outside  the  sub-domain  taken  for  the 
polar  integral.  But,  for  open  surfaces,  it  is  inevitable  to  have  non- vanishing 
integrand  if  the  target  point  is  near  the  open  edge.  We  resolved  this  problem  by 
re-sampling  data  with  Chebyshev  spacing  from  the  uniformly  distributed  data 
for  each  angle,  if  the  radial  line  cuts  the  open  edge. 

-  The  radial  integral  as  a  function  of  angle  exhibits  singular  behavior  as  the  tar¬ 
get  point  approaches  to  the  open  edge.  To  resolve  this  singularity >  we  applied 
polynomial  sigmoidal  transformation  for  the  angle  variable. 

•  At  first  our  solver  experienced  an  instability  during  the  solve  of  linear  system  arising 
from  TV-  We  suspect  that  the  instability  results  from  the  following  facts:  (1)  the 
second  integral  in  T\y  is  a  pseudo-differential  operator  of  order  L  Repeated  application 
of  such  an  operator  can  generate  growing  oscillating  errors.  (2)  Unlike  scalar  problems, 
our  operator  involves  both  1/Vd  and  Vd  scaling  in  a  single  equation.  A  small  error  in 
au  and  av  might  result  in  inadequate  density  scaling.  We  resolved  this  instability  by 
introducing  a  regularization  of  the  form  T\ySwm  ■  where  S\ym  is  a  single  layer  potential 

with  scaling  matrix  W*  —  ^  1  ^ 


The  figure  shows  the  result  of  our  first  successful  implementation  of  this  method.  The 
method  works  correctly  and  our  regularization  successfully  stabilized  the  linear  solver;  we 
confirmed  that  our  formulation  is  correct  and  appropriate.  In  the  next  phase  of  our  re¬ 
search  in  this  context,  we  will  focus  on  (1)  Implementation  of  well-engineered  version  of  the 
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Figure  6:  Detached  and  re-attached  pylon,  front  view.  Left:  Original  low  quality  surface 
rendering.  Right:  processed  surface  by  means  of  the  continuation  method  surface  represen¬ 
tation  algorithm. 


code,  (2)  Enhancement  of  the  preconditioner  to  achieve  further  reductions  in  the  number 
of  iterations  required  by  the  linear  algebra  (GMRES)  iterative  solver. 


2  Enhancement  and  Extension  of  Geometry  Representation 
Software 

Significant  enhancements  were  introduced  to  the  geometry  representation  capability,  en¬ 
abling  treatment  of  heavily  damaged  structures.  The  new  capabilities  are  demonstrated  in 
Figure  1  and  Figures  5  through  8.  Clearly,  the  new  software  can  repair  surfaces  containing 
very  challenging  features,  including  very  coarse,  low  quality  discretizations,  detachment, 
lack  of  water-tight  ness,  etc. 


3  Simplified  methodologies  for  treatment  of  Corners  and  Edges 

Our  strategy  for  simplification  of  problems  concerning  corners  and  edges  is  two-fold:  on  one 
hand  we  developed  1)  Improved,  highly  accurate  solvers  for  the  accurate  solution  of  such 
problems  even  in  cases  in  which  the  surface  currents  and  physical  fields  tend  to  infinity,  and, 
we  then  produced  2)  A  simpler,  more  easily  implementable  approximate  approach— whose 
accuracy  was  tested  for  the  Dirichlet  boundary  value  problem  in  Section  3.2  by  comparison 
with  a  rigorous,  previously  existing  highly  accurate  method.  We  describe  these  two  method¬ 
ologies,  in  turn,  in  the  next  two  subsections.  The  advantage  of  the  simplified  strategy  2) 
is  that  it  does  not  require  complicated  special  treatment  of  the  corner  (like  that  presented 
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Figure  7:  Detached  and  re-attached  pylon,  side  view;  see  also  Figure  6.  Left:  Original  low 
quality  surface  rendering.  Right:  processed  surface  by  means  of  the  continuation  method 
surface  representation  algorithm. 

in  Section  3.1),  or  the  much  more  complex  thxee-d  conical  points),  for  which  solution  of 
eigenvalue  problems  is  necessary  to  evaluate  the  blow-up  exponents.  In  subsequent  work  we 
will  carry  out  the  similar  comparisons  for  the  Neumann  boundary  value  problem  (for  which 
the  electric  field  blows  up,  and  for  which  the  exact  solution  will  be  obtained  by  means  of 
the  numerical  method  mentioned  in  point  1)  above),  at  which  point  an  extension  to  the  full 
three-d  Maxwell  problem  should  prove  rather  direct. 


3.1  Corners  and  Edges  and  infinite  solutions:  rigorous  formulation 

We  have  introduced  a  new,  highly  improved,  highly  accurate  algorithm,  based  on  integral 
equation  formulations,  for  the  treatment  of  Partial  Differential  Equations  in  closed  two- 
dimensional  domains  with  non-smooth  boundaries;  we  focus  in  cases  in  which  the  integral- 
equation  solutions  as  well  as  physically  meaningful  quantities  (electric/magnetic  fields)  tend 
to  infinity  at  singular  boundary  points  (corners);  see  e.g.  Figure  9.  Our  numerical  results 
demonstrate  excellent  convergence  as  discretizations  are  refined,  even  around  singular  points 
at  which  solutions  tend  to  infinity.  We  believe  this  is  the  first  contribution  in  which  infinite 
solutions  at  corners  and  edges  arc  captured  accurately.  The  main  elements  of  our  approach 
are  described  in  what  follows. 


Asymptotic  expansions  and  cancellation  of  infinities.  The  essence  of  our  method 
resides  in  explicit  cancellation  of  infinities  that  arise  as  asymptotic  expansions  of  the  form 

<p{r,  9)  -  ar~$  +  o(l),  ,J 3  >  0,  (6) 
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Figure  8:  Full  aircraft*  Left:  Original  low  quality  surface  rendering*  Right:  processed 
surface  by  means  of  the  continuation  method  surface  representation  algorithm, 

(where  r  is  the  polar  distance  around  an  edge  x  =  {0}<-*{r  =  0})  are  substituted  in  an 
integral  equation  of  the  form 

f  ^^(V)  dS{y)  =  s(x)  for  x  G  T  \  {0}  ,  f  <j>(y)  dS{y)  -  0  ;  (7) 

the  (smoother)  difference  between  the  exact  solution  and  the  asymtpotics  (6)  can  then  be 
obtained  by  means  of  previously  introduced  integral  equation  solvers,  In  detail,  note  the 
indeterminate  limits  that  arise  as  t  —>  0  and  t  —*  T — since  (p  —+  oo  at  the  edge,  and  thus  so 
must  the  integral,  in  order  to  agree  with  the  finite  limit  of  the  right  hand  side.  To  eliminate 
the  numerical  cancellations  and  accuracy  loss  that  would  otherwise  occur  for  t  close  to  0 
and  T  we  resolve  in  closed  form  this  cancellation  of  infinities,  as  explained  in  what  follows. 

Remark  3.1  .4  numerical  illustration  study  demonstrating  the  need  for  a  special  numerical 
treatment  of  this  problem  has  been  undertaken:  we  have  concluded  in  no  uncertain  terms 
that  no  accuracy  may  be  expected  in  the  integral  equation  solution,  in  a  neighborhood  of  the 
corner  or  edge *  unless  this  issue  is  satisfactorily  addressed. 

To  resolve  the  indeterminate  limit  mentioned  above  we  first  re-express  the  integral  equa¬ 
tion  solutions  to  explicitly  account  for  the  asymptotics  (6):  we  write 

y(t)  =  pi(t)  +/i2(t)  >  MO  =  ,  q  =  (  1  _V  ’  J  <  2  '  ^ 

V  Q  * 
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Figure  9;  Obtuse  and  acute  corners, 
where  ip  is  a  smooth  real-valued  function  for  0  <  t  <  T  satisfying 


'1  for  t  e  [0.  fi]  U  [S,T]  when  0  <  a  <  1,  and 


ip{t)  =  < 


—  1  for  t  €  (0,  fi] 
1  for  t  €  [S,  T] 


when  1  <  a  <  2, 


and  where  /X2  is  a  Holder  continuous  function;  the  unknowns  in  our  problem  thus  become 
the  parameter  a  €  R  and  the  function  ^2-  The  integral  equation  can  now  be  expressed  in 
the  form 


-  f  K(t,s)p2(s)ds  +  afi(t)  =  f(t)  ,  [  ja;'(s)|^2(s)ds  +  az\  =  0.  (9) 

^  J  o  Jo 

where 

AW  =  -  [  K(t,s)fj,1(s)d$  and  zi  =  [  |aj'(s)|jUi(s)  ds.  (10) 

*  Ju  J  0 

Clearly,  the  infinity  cancellation  mentioned  in  Remark  3.1  is  now  encapsulated  in  the  func¬ 
tion  fi  in  (10) — since  the  function  \i2  is  bounded,  and,  thus,  so  are  all  other  terms  in 
equation  (9). 


Stable  evaluation  of  indeterminate  limits.  The  difficulties  arising  from  cancellations 
of  infinities  can  be  addressed  by  resolving  analytically  the  indeterminate  limits  lim^Q  fi  (0 
and  lim^x /i(t).  To  do  this,  we  first  introduce  approximations  L\  and  L2  of  the  kernel  K 
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around  the  endpoints  0  and  T\ 


(ta'(O)  -  (s  -  7V(r))  ■  n(0)  |q-'(T)|  sin(«7r)  |(.s  -  T)sf{T)\  |s'(r)| 

U  ’  }  2Tv\tx'(0)-{s~T)x'(T)\2  2Tr\tx’(0)-{s-T)x'(T)\2  [  1 


r  ( f  A  =  ((*  ~  T)x\T)  ~  ax>(0))  •  n(T)  =  sin(Q7r)  |.sr,(Q)|  |a;'(0)| 

21  ’  1  2tt| (t  -  T)x'{T)~  sx'(0)|2  '  2tt| (t-T)x'(T)  -  sx'(0)|2  ' 


and  the  associated  quantities 


(13) 


(14) 


Remark  3.2  IfTi  and  r2  are  straight  segments  parametrized  by  linear  functions  with  speeds 
xf(Q)  arid  x^T),  then  L\(t,  s )  —  K(t ,  s)  for  t  near  0  and  s  near  T,  and  L2(t ,  s)  =  K(t ,  s)  for 
t  nearT  and  s  near  0.  The  quantities  a\  and  &2 ,  which  themselves  give  rise  to  cancellation 
of  infinities  as  t  -4  0  and  t  —*  T ,  capture  the  essence  of  the  cancellations  associated  with 
the  quantity  fi(t). 

The  key  to  our  algorithm  lies  in  recognizing  that  <i\  and  a 2  can  be  evaluated  as  rapidly 
converging  series  expansions,  thereby  providing,  via  the  simple  additional  manipulations 
described  in  Section  3-1,  an  efficient  and  numerically  stable  means  for  evaluation  of  the 
quantity  f\(t)  in  and  mound  t  —  0  and  t  —  T. 

Series  expansions  for  a\  and  02-  Introducing  the  notations 


U  =  cr-s)|x'(T)| 


V2  =  ii|x'(0)| 


together  with  the  parameters  0  <  S  <  R  <  T.  which  are  selected  in  such  a  way  that 
B\ {t)  <  1  for  t  <  S,  and  B2(£)  <  1  for  t  >  R ,  we  have  the  following  series  expansions  for 
the  functions  and  cj2: 


for  t  G  (0,3}  (15) 


for  t  €  [/i,T).  (16) 


10 


In  particular.  <J\  and  02  are  analytic  functions  on  their  domains  of  definition.  We  can.  in 
fact,  establish  a  somewhat  more  general  result  for  the  functions 


vf(t,  A)  =  ±^T  , 

cr±(t,  X)  -  — — Tf  ^  ±  jf  L2(i,s)jsx'(0)jAds  . 


Lemma  3*3  Let  —1  <  \  <0  be  given ,  For  t  £  (0, 5],  we  have 


*f(M)  =  (i  ±  si°(^:((Lr)t,,r))  ^  T I  £  (IT) 


cma/ogtms  e253a.ns.i0n  ftdds  /or  m  £he  domain  t  €  [^,T), 

Remark  3*4  In  light  of  the  facts  that }  ai(t)  =  q)  and  <12 (0  —  g)  /or  0  < 

q  <  L  ru/ide  o" i  (£)  =  — er+(t,  — g)  and  ct2 (^ )  —  (t,  — g)?  /or  1  <  a  <  2,  the  reiataons  (15) 
and  (16)  /oiio^  dzrec%  /rom  the  previous  lemma  -since,  for  A  —  —q  (with  q  given  by  (8) ), 
the  coefficient  0/ |£a/(0)|A  in  (17)  vanishes. 

Evaluation  of  the  remainders  f\(t)  —  (t)  and  f\{t)  —  02  (i)*  Having  computed  crj(t) 

and  craft),  the  evaluation  of  /i  (t)  reduces  to  evaluation  of  certain  singular  integrals  (that  can 
be  treated  numerically  by  previously  existing  methods  and  generalizations  thereof)  together 
with  some  simple  algebraic  manipulations.  Indeed,  to  evaluate  /i(t)  we  define 


E^s)  =  Li{t,s)\xf{T)\-*  -  K{t,s) 


x(s) 


s-T 


~<1 


,(1)  =  #)ffitMor=il(0) 


x(t) 


-q  \  t-q 

-wm~q  V 


E2{t,s)=if>{0)  L2(t,s)\x'(0)\-<>  -  K(t,s) 


x(s) 


and 


m(t)  = 


x(t) 


t-T 


(18) 

(19) 

(20) 

(21) 


and  we  employ  three  different  formulae,  one  for  each  of  the  cases  t  G  [0,5],  t  €  (S,  R)  and 
t  G  \R:T]. 
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For  t  e  [0,  S]  we  use 


/l(t)'-=^i(*)  +  ffi(0-  fS  K{t,s)fiib)ds+  f  E\{t,  s){T  —  s)~q  ds  .  (22) 

Jo  Js 


For  t  6  (5,  R),  we  use 


A  C0  =  ^p-j( 


(23) 


Finally,  for  t  G  [R,T],  we  use 


/l(0  =%(*)  + ff3 (^3  “  f  K{t,s)fjii(s)ds+  f  E2(t,s)s  q ds  .  (24) 

Jh  Jo 


It  is  easy  to  check  from  equation  (19)  that  /71  (t)  is  bounded  throughout  [0»T].  The 

integral  in  (22)  is  clearly  bounded  for  t  G  [0,5],  and  so  is  the  quantity  fJ1(t).  It  follows 

that  the  integral  in  (22)  bounded  for  t  G  [0,5],  and  our  procedure  has  thus  adequately 
resolved  the  indeterminate  limit  for  equation  (7)  as  t  — >  0;  similarly  the  approach  above 
resolves  the  indeterminate  limit  as  t  — *  T. 

The  quantities  7}\  (£),  Tfy{t)  can  be  computed  directly  using  the  expressions  (19)  and  (21) 
for  t  sufficiently  far  from  0  and  T.  In  order  to  avoid  cancellation  errors,  in  turn,  for  t  close  to 
0  our  algorithm  for  evaluation  of  771  (t)  uses  a  few  terms  of  the  Taylor  expansion  of  the  term 
in  parenthesis  in  the  right-hand  expression  of  equation  (19);  similarly,  a  Taylor  expansion 
is  used  to  evaluate  772  (t)  for  t  close  to  T. 

This  algorithm  has  given  rise  to  solutions  with  full  double  precision  accuracies  in  near 
fields,  far  fields  and  surface  currents,  for  highly  challenging  domains  with  corners  such  as 
those  depicted  in  Figure  9.  We  believe  this  is  the  first  algorithm  capable  of  resolving  the 
challenging  Neumann  problem  with  high  accuracy  around  corners  and  edges;  a  solution  for 
the  significantly  simpler  Dirich  let  problem,  in  turn,  had  been  provided  by  R.  Kress  and 
other  authors,  and  has  been  available  for  a.  number  of  years. 

3.2  Corners  and  Edges  and  infinite  solutions:  simplified  formulation 

While  the  approach  introduced  in  Section  3,1  is  perfectly  accurate  and  reliable,  its  im¬ 
plementation  is  somewhat  more  complicated  than  might  be  desirable:  we  thus  pursued  an 
additional,  alternative  avenue,  in  which  the  geometry  is  smoothed  prior  to  solution  of  the 
integral  equation  problem.  Our  strategy  thus  is  to  use  highly  accurate  solutions  for  singular 
geometries,  such  as  those  developed  in  the  previous  section,  to  answer  the  following  ques¬ 
tion:  what  smoothings,  if  any,  can  be  used  while  giving  accurate  near  and  far  field  values? 
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s  =  1 

tto 

1! 

W 

C* 

II 

-A 

s  =  15 

d-  0.01 

d  ~  0.005 

d  ~  0.0025 

d  ~  0.001 

Table  1:  Maximum  distance  (equal  to  the  distance  between  the  actual  corner  and  the 
smooth  approximating  curve)  for  various  values  of  the  sharpness  paramteter  s. 


How  close  to  the  singularity  can  any  accuracy  be  expected  from  the  smoothed  edge  solver? 
We  have  thus  far  answered  (satisfactorily)  these  questions  in  the  case  in  which  the  solution 
does  not  tend  to  infinity  at  the  corner;  a  study  the  complementary  case  is  being  initiated 
as  of  this  writing  and  should  be  completed  shortly. 


Figure  10:  Sharp  (blue)  and  smoothed  (red)  corner,  sharpness  level  s  =  1  (smoothest 
approximation;  other  smoothings  we  used,  see  Tables  1  and  2,  give  rise  to  much  sharper 
approximations  of  the  corner.  Left:  full  curve.  Right:  corner  closeup. 

The  result  of  our  study  on  the  effect  of  smoothing  is  presented  in  Figure  10  and  Ta- 
hie  2;  a  visualization  of  the  surface  representation  methodologies  we  developed  to  produce 
(automatically)  smoothing  of  corners  and  edges  is  shown  in  Figure  11, 

In  Figure  10  we  depict,  in  blue,  a  scatterer  containing  a  sharp  corner  together  with 
a  smooth  approximation  of  it,  shown  in  red.  The  smooth  curve  depicted  in  this  figure  is 
the  smoothest  (i*e*,  the  one  with  the  most  rounded  corner)  in  a  series  of  approximating 
smooth  curves  (determined  by  the  value  of  a  “sharpness  parameter"  s )  considered  in  this 
section:  see  Table  L  For  example,  the  largest  distance  d  between  the  exact,  sharp-corner 
curve  and  this  s  —  1  approximation  is  d  —  0.01;  the  corner  distances  for  the  various  smooth 
approximations  we  considered  are  displayed  in  Table  1. 

In  Table  2,  in  turn,  we  consider  solutions  of  a  Dirichlet  problem,  for  which  solutions 
remain  bounded,  and  we  present  the  errors  that  arise,  at  various  points  in  space,  as  the 
electric  field  for  the  sharp-corner  scatterer  is  approximated  by  the  smooth  corner  scatterer. 
Clearly  the  approximations  are  excellent  in  the  far  field,  and  they  provide  quantitatively 
accurate  approximations  even  and  and  very  close  to  the  comer  point  itself *  The  sharp  corner 
solution  was  computed  by  means  of  a  highly  accurate  Dirichlet  solver;  we  will  soon  present  a 
similar  approximation  procedure  for  the  Neumann  problem  (for  which  the  electric  field  blows 
up,  and  for  which  the  exact  solution  will  be  obtained  by  means  of  the  numerical  method 
described  in  Section  3*1),  at  which  point  an  extension  to  the  full  three-d  Maxwell  problem 
should  prove  rather  direct.  The  main  advantage  of  the  corner  smoothing  approach  is  that, 
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Sharpness  level 

Max  Error 

(0,0) 

(-0.1,0) 

(-1,0) 

(-10,0) 

s  -  1 

1.8  xlO-2 

3.0  x  10-3 

1.2  x  10:J 

9.2  xlO-4 

s  =  3 

i.i  xio-* 

1.2  x  10-3 

3.9  x  10-4 

3.0  xlO"4 

s  =  7 

7.1  xl0~3 

4.3  x  10-4 

1.4  x  10“ 4 

1.1  xlO-4 

s  =  15 

5.7  xlO-3 

3.2  x  10-4 

4.9  xlO"5 

3.5  x  10_a 

Table  2:  Errors,  at  various  points  in  space,  in  the  approximation  of  the  true  field  (that  is. 
the  field  scattered  by  the  scatterer  with  the  perfectly  sharp  corner)  by  the  field  scattered 
by  the  smoothed  scatterer.  The  corner  itself  is  at  the  point  (0,0).  In  all  cases  a  total  of 
128  discretization  points  were  used ,  both  in  the  solution  for  the  perfectly  sharp  corner  and 
for  the  smoothed  geometries.  The  smoothest,  less  accurate  approximation  {sharpness  level 
s  =  1)  is  depicted  in  Figure  10.  Note  that  convergence  is  observed,  evan  at  the  corner  itself; 
very  high  accuracies  result  in  the  all  important  far  field  values. 


clearly,  no  complicated  special  treatment  of  the  corner  (like  that  presented  in  Section  3.1), 
or  the  much  more  complex  three-d  conical  points)  is  needed. 


Figure  11:  Demonstration  of  the  surface-smoothing  capability.  Left:  full  smoothed  surface. 
Right:  corner  doseup. 


4  Integration  of  surface  representation  and  solver  capabili¬ 
ties:  application  to  realistic  structures. 

Additional  work  over  the  last  performance  period  focused  on  integration  of  the  surface 
representation,  solver  and  GUI  capabilities.  The  results  of  this  (still  ongoing)  effort,  are 
shown  in  Figure  1  Right,  As  demonstrated  in  this  figure,  the  geometry  representation 
algorithm  for  smooth  surfaces  as  been  integrated  with  the  salvers;  after  completion  of  the 
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solvers  based  on  edge-smoothing  (as  described  in  Section  3.2) ,  complete  software  integration 
will  be  achieved  and  a  capability  for  solutions  of  unrestricted  geometries  from  CAD  will 
be  available.  To  accomplish  this  goal,  in  particular,  it  will  be  necessary  to  produce  a 
methodology  enabling  easy  smoothing  of  corners  and  edges,  with  arbitrarily  sharpness  levels, 
as  described  in  Section  3.2.  This  effort  is  already  in  progress;  preliminary  results  of  the 
resulting  software  are  presented  in  Figure  11. 

5  Curved  wire  antenna  solvers 

We  developed  an  extension  of  the  methodology  previously  put  forward  in  collaboration 
with  M,  Haslam  for  straight  wire  antennas  (“Regularity  theory  and  super- algebraic  solvers 
for  wire  antenna  problems”,  O,  Bruno  and  M,  C,  Haslain;  SIAM  Jour,  Sci.  Comp.  29 
1375-1402,  (2007))  that  applies  to  arbitrary  curved  wire  antennas.  Thus  far  the  code  has 
not  been  optimized,  and  is  restricted  to  closed  curved  wires  driven  by  an  incident  field;  an 
extension  to  arbitrary  curved  wires  and  to  driving  sources  on  the  wire  itself  is  the  subject 
of  our  continuiung  work  in  these  regards,  (We  mention  that  the  fundamental  problem  of 
incorporating  driving  sources  with  high-order  accuracy  was  dealt  with,  for  straight  wires, 
as  part  of  earlier  work,)  The  present  curved  wire  code  is  based  on  use  of  the  electric  field 
integral  equation  with  a  thin- wire  approximation — according  to  which  the  electrical  current 
J  is  constant  around  the  circular  cross-section  that  results  as  the  wire  is  cut  by  an  orthogonal 
plane.  This  approximation  is  commonly  used  in  the  literature;  for  example  we  mention  the 
contributions  provided  by  Wu  (1962)  for  a  circular  geometry,  as  well  as  the  recent  article 
by  Champagne  and  Wilton  (2006).  (The  latter  paper  incorporates  high-order  elements; 
however,  since  the  logarithmic  singularity  is  treated  using  a  tangent  line  approximation  the 
approximation  is  reduced  to  lowest  order,)  We  believe  our  algorithm  is  the  first  one  to  deliver 
high-order  accuracy  for  the  curved- wire  problem.  Our  preliminary  results  are  as  follows; 
for  a  circular  wire  with  circle  radius  given  by  kb  =  5,  and  wire  thickness  ka  —  0.01  (i,e.,  a 
challenging  very-thin  wire  problem),  our  code  produces  solutions  with  one  percent  error  in 
a  small  fraction  of  a  second  in  a  single  processor,  while  a  couple  of  seconds  calculation  give 
very  many  digits  of  accuracy. 

Personnel  Associated  with  the  Research  Effort 

A,  Anand.  O.  Bruno,  T,  Elling,  M.  Haslam,  J.  Qvall,  F.  Reitich  and  C.  Turc. 


Publications 

*  O.  P.  Bruno,  Jeffrey  S,  Ovall,  and  Catalin  Turc  i:A  High-Order  Integral  Algorithm 
for  Highly  Singular  PDE  Solutions  in  Lipschitz  Domains”;  In  preparation. 

#  O.  Bruno  and  Andy  Monro  J.,  “Spectrally  accurate  windowing  for  rough-surface  scat¬ 
tering  problems”;  In  preparation. 

•  O.  Bruno  and  Stephane  Lintner,  “A  highly  accurate  integral  method  for  open-surface 
scattering  requiring  small  numbers  of  GMRES  iterations”:  In  preparation. 

*  O.  Bruno  and  Stephane  Lintner.  uWell-posed  Integral  Equations  for  Acoustic  Prob¬ 
lems  on  Open  Surfaces” ,  In  preparation. 
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•  0,  P.  Bruno  and  M.  Haslam,  “High-Order  Solution  of  the  Scattering  Problem  for 
One- Dimensional  Perfectly- Conducting  Periodic  Surfaces”;  Submitted, 

•  O.  Bruno,  T.  Elling,  R.  PafFenroth  and  C.  Turc,  “Electromagnetic  integral  equations 
requiring  small  numbers  of  GMRES  iterations”.  Submitted, 

•  J.  Chaubellj  C.  O.  Ao  and  O.  P.  Bruno,  “Evaluation  of  Propagation  of  GPS  Signals 
through  the  Atmosphere  via  High-Frequency  Localization  and  Rytov’s  Approxima¬ 
tion”,  To  appear  in  Radio  Science 

•  “Superlens-cloakmg  of  small  dielectric  bodies  in  the  quasistatic  regime”  O .  Bruno  and 
Stephane  Lintner,  Journal  of  Applied  Physics  102,  124502  (2007) 

•  “Accurate,  high-order  representation  of  complex  three-dimensional  surfaces  via  Fourier- 
Continuation  analysis” ,  O.  Bruno,  Y.  Hari  and  M,  Pohlman;  Journal  of  Computational 
Physics  227  (2007)  1094-1125 

•  “Regularity  theory  and  super-algebraic  solvers  for  wire  antenna  problems”,  O.  Bruno 
and  M.  C.  Haslam;  SIAM  Jour.  ScL  Comp.  29  1375-1402,  (2007) 

•  “Time  Stepping  Via  One-Dimensional  Fade  Approximation”,  D,  Amundsen  and  O. 
Bruno;  J.  Sci.  Comput.  30  83-115.  DOI  10. 1007/sl091 5-005-9021-4,  (2007) 

•  “An  0(1)  Integration  Scheme  for  Three-Dimensional  Surface  Scattering  Problems”, 
O.  Bruno,  and  C.  Geuzaine;  Journal  of  Computational  and  Applied  Mathematics  204 
463-476  (2007). 

Interactions /Transitions 


Participation  in  a  wide  range  of  scientific  meetings,  interaction  with  university  colleagues 
and  seminars  at  universities,  as  well  as  seminars  and  consulting  work  for  industrial  parties 
are  a  regular  part  of  our  research  efforts.  Recent  and  forthcoming  plenary  lectures  are 
listed  in  the  section  “Honors  and  awards”  below.  Our  most  recent  industry  interactions 
concerns  work  on  optics  problems  arising  from  “New  planet  discovery”,  in  collaboration  with 
Northrop  Grummann  researchers  Amy  Lo  and  Tiffany  Glassman,  and  collaborations  with 
a  JPL  (NASA)  group  (J.  Chaubell,  B.  Pollard  and  E.  Rodriguez)  concerning  development 
of  accurate  scattering  solvers  applicable  to  large  reflectarray  geometries. 

Honors/ A  wards 


Oscar  Bruno  was  plenary  speaker  at  the  following  conferences. 

*  “TBA”.  plenary  speaker  at  “ICOSAHOM  2009”  (International  Conference  on  Spectral 
and  High  Order  Methods)  http://svein.halvorsen.ee/icosahoiii/frog/71iome/. 

#  “Integral  equations  in  regular  and  singular  domains”,  in  the  Gberwolfach  meeting 
on  “Analysis  of  boundary  element  methods”;  Martin  Cotabel  and  Ernst  P.  Stephan, 
Organizers,  April  13th  -  19th,  2008. 
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•  “Fast  High-Order  High-Frequency  Solvers  in  Computational  Acoustics  and  Electro¬ 
magnetics”  ,  in  Seventh  Mississippi  State  -  UAB  Conference  on  Differential  Equations 
&  Computational  Simulations  November  1-3,  2007  Doubletree  Hotel  Birmingham, 
AL,  USA.http:  //www.msstate.  edu/dept/math/events/de,  conf /de2007/ 

•  “Fast  High- Order  High-Frequency  Solvers  in  Computational  Acoustics  and  Electro¬ 
magnetics”  ,  in  Seventh  Mississippi  State  -  UAB  Conference  on  Differential  Equations 
&  Computational  Simulations  November  1-3,  2007  Doubletree  Hotel  Birmingham, 
AL,  US  A. http://iifww.ms  st  ate .  edu/dept /math/event  s/de .  conf /de2007/ 

•  Organizer,  in  collaboration  with  R.  Kress,  of  the  workshop  “High-order  methods  for 
computational  wave  propagation  and  scattering”  at  the  American  Institute  of  Math¬ 
ematics  (Palo  Alto,  CA). 

http : / /www . aimath , org/ARCC/workshops/wavescattering . html 

•  “High-order  scattering  solvers  and  surface  representation  algorithms”  invited  Talk  at 
ICIAM’s  Industry  Days  “Computational  Electromagnetics”.  July  16-20  2007,  07. 
http : / / www , ic i am07  * ch/embedded_me  e t i ngs / 1 ndus  t ry Day s 

•  “Fast  High- Order  High-Frequency  Solvers  in  Computational  Acoustics  and  Electro¬ 
magnetics”;  leading  lecture  at  the  Isaac  Newton  Institute  Workshop  Effective  Com¬ 
putational  Methods  for  Highly  Oscillatory  Problems:  The  Interplay  between  Mathe¬ 
matical  Theory  and  Applications  2  July  to  6  July  2007;  A  week-long  workshop  at  the 
Isaac  Newton  Institute  20  Clarkson  Road,  Cambridge  CB3  OEH,  UK. 

•  Participant  at  Cambridge  University's  Newton  Institute  as  part  of  the  program  on 
Highly  Oscillatory  Problems;  June/ July  2007. 

•  “Efficient  evaluation  of  high-frequency  propagation  and  scattering,  with  application  to 
propagation  in  non-spherically-symmetric  atmospheres”,  in  Oscillatory  integrals  and 
integral  equations  in  high  frequency  scattering  and  wave  propagation”,  June  19,  2007. 
A  one-day  workshop  at  the  Isaac  Newton  Institute  20  Clarkson  Road,  Cambridge  CB3 
OEH,  UK. 

•  “Computational  electromagnetics  and  acoustics:  high-order  solvers,  high-frequency 
configurations,  high-order  surface  representations”;  45  min,  leading  presentation  at 
the  Oberwolfach  meeting  on  Computational  Electromagnetism  and  Acoustics,  Febru¬ 
ary  5  to  9th,  2007. 

•  “Fast  and  accurate  electromagnetics  simulations”  (at  ETOPIM  7:  The  7th  Interna¬ 
tional  Conference  on  the  Electrical,  Transport  and  Optical  Properties  of  Inhomoge¬ 
neous  Media.  Sydney,  Australia.)  July  9-13,  2006 
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